Frustrated Heisenberg antiferromagnet on a pyrochlore lattice 
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We investigate quantum phase transitions for the s = 1/2 antiferromagnetic Heisenberg model on 
a pyrochlore lattice. By means of a series expansion starting from isolated tetrahedra, the ground- 
state phase diagram is determined. When the ratio of the two competing exchange couplings is 
varied, the first-order (second- order) quantum phase transition occurs between the two spin gap 
phases (the spin-gap and the antiferromagnetic phases). We also discuss some properties expected 
for the s = 1 pyrochlore spin system. 



Geometrically frustrated magnetic materials have been 
the subject of considerable interest recently. A typical ex- 
ample is a class of transition-metal oxides with pyrochlore 
structure. In particular, it was reported that the specific- 
heat coefficientJs exceptionally large for a metallic com- 
pound LiV 2 04.uIt was recently pointed out that the frus- 
tration caused by a tetrahedral network of vanadium ions 
may be important to understand the heavy-fermion be- 
havior in this compoundp stimulating further intensive 
studies on the related systems with pyrochlore structure. 
Such frustration effects should be much morejjuximinent 
for quantum spin systems. Theoretical studiesautj for the 
s = 1/2 quantum spin model on a pyrochlore lattice were 
first done by Harris et al. p who pointed out the possibil- 
ity of the dimerized ground state by exploiting a field the- 
oretic approach. Canals and Lacroixel clarified that the 
ground state of the model is a spin-liquid state with the 
spin gap. They found that the neutron diffraction cross 
sectionO observed ip. Y(Sc)Mn 2 is in fairy good agreement 
with their resultsia Isoda and Morija however, used a 
bond-operator approach to suggest that the ground state 
may be described by a RVB-like tetrahedral (plaquette) 
singlet state, which is different from the dimer-singlet 
state known so far. Furthermore, the possibility of the. 
"topological spin glass" was pointed out in Y 2 Mo 2 C-7,lI 
making this issue more attractive and challenging. 

In this paper, we investigate the s — 1/2 quantum spin 
model on a pyrochlore lattice with competing antiferro- 
magnetic interactions shown in Fig. [l] (a). Our system 
may describe some pyrochlore-lattice compounds such as 
Y(Sc)Mn 2 la as well as GeCu 2 4 el found recently. We 
study the ground-state phase diagram and clarify the role 
of the geometrical frustration by studying quantum phase 
transitions by means of series expansion techniques^ We 
also discuss the s = 1 case briefly. 

Let us first consider the s = 1/2 spin model on a 
pyrochlore lattice, which is described by the following 
Hamiltonian 



FIG. 1: (a) Frustrated antiferromagnetic spin model on a 
pyrochlore lattice. Bold and thin-solid lines represent the 
exchange couplings J and J', respectively, (b) Initial config- 
uration for the series expansion: Broken lines represent the 
perturbed bonds, A J and A J' (see text). 
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where (i, j) denotes a pair of two adjacent sites connected 
by the thick bond in Fig. [l] (a), whereas < i, j > is that 
for the thin bond. Both of the exchange couplings J 
and J' are assumed to be antiferromagnetic. For the 



compound Y(Sc)Mn 2 , we may .take J/J' ~ 1,0 while for 
GeCu 2 4 , it is to be J/J' ~ 6.1 

Before proceeding with the analysis, we note that, in 
some limiting cases, the spin system is reduced to sim- 
ple models whose nature can be easily understood. For 
J' = 0, it is equivalent to the s — 1/2 massless Heisen- 
berg spin chain, whereas an antiferromagnetically or- 
dered state is stabilized for J = since the model has 
a three-dimensional structure without frustration in this 
case. To determine the phase diagram, we study what 
kind of quantum phase transition occurs when the com- 
peting interactions are varied— For this purpose, we use 
the series expansion methodjj which has an advantage 
to deal with frustrated spin systems in higher dimen- 
sions. In fact, it was successfullyi-applied to frustrated, 
spin systems such as Ji — J*,model,ll3 plaquette system,£ll 
orthogonal-dimer system.E3 To apply the series expan- 
sion method to the pyrochlore-lattice system, we first 
divide the original Hamiltonian eq. (Q) into two parts 
as H = Hq + XHi by introducing an auxiliary param- 
eter A, where Hq(Hi) represents the unperturbed (per- 
turbed) Hamiltonian. Note that the system is reduced to 
the original model for A = 1. We here choose a tetrahe- 
dron [Composed of four spins as a starting configuration 
(i/ ),&El and then connect each tetrahedron via tetrahe- 
dral bonds labeled by AJ and AJ' [see Fig. | (b)]. The 
Hamiltonian h for an isolated tetrahedron in Hq is given 
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by 

h = J (Si • S 2 + S 3 • S 4 ) + J' (Si + S 2 ) ■ (S 3 + S 4 )(2) 

The energy eigenvalues, E, of the tetrahedron for a given 
j = J/J' are listed in Table B, where 6*12 (5 I 34) represents 
the combined spin Si + S 2 (S3 + S4), and Stotai is the 
total spin. It is seen in this table that for < j < 1, 



TABLE I: Eigenvalues of the s = 1/2 spin system on an iso- 
lated tetrahedron. 
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the isolated tetrahedron has a plaquette-singlet ground 
state with S12 = S34 = 1 and Stotai = 0. On the other 
hand, for j > 1 we have the dimer ground state with 
S12 = S34 = Stotai = 0. The phases specified by these 
singlets are referred to as the plaquette and dimer phases, 
respectively. 

Keeping the above properties in mind, we now discuss 
how the plaquette and the dimer states compete with 
each other when the inter-tetrahedron couplings A J and 
AJ' are introduced. We expand the ground state energy 
up to the sixth order in A for several values of j. We 
show the obtained energy in Fig. ||, for jaihich the first- 
order inhomogeneous differential methodtij is applied to 
the bare series. As mentioned above, for A — > 0, the first- 



-0.3 
-0.4 
Eg/NJ' -0.5 
0.6 



0.7 



k=0.0 ,.. -><;v 












^X=1.0 







0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

j 



FIG. 2: Ground state energy for the dimer (right-side) and 
the plaquette (left-side) phases for various values of A. The 
energy for A = 0.0,0.5 and 1.0 is shown as the dot-dashed, 
dashed and solid lines, respectively. 

order quantum phase transition occurs between the pla- 
quette and the dimer phase at the critical point j c = 1. 
It is seen in the Fig. || that the critical value for the 
phase transition is little changed even if we increase A. 
In fact, the energy up to the second order in A is same 
for_both states near j = 1, as pointed out by Harris et 
al.u For small A, the first-order transition point is given 



by jc ~ 1 — 0.021A 3 . Remarkably enough, the first-order 
transition point is estimated as j c ~ 1.0 even for A = 1, 
where our generalized model is reduced to the original 
one. Therefore, we arrive at a quite interesting conclu- 
sion that the homogeneous spin system (A = j = 1) 
with pyrochlore structure is located quite closely to the 
phase boundary of the first-order quantum phase transi- 
tion, although it is difficult to definitely conclude which 
phase the ground state really belongs to within our ac- 
curacy. This fact clarifies the reason why Harris et al. 
and Isoda et al. had different conclusions on the nature 
of the ground state for the same j = 1 model, where the 
former (latter) claimed that the ground state is a dimer 
singlet (plaquette singlet). As mentioned above, the en- 
ergy for two phases is very close to each other, so that 
the mean-field type treatment may not correctly specify 
the ground state. Furthermore, there even remains the 
possibility that the system is just on the boundary, and 
thus the ground state could be degenerate at j = 1. In 
any case, it is instructive to notice that unusual dual- 
properties reflecting both natures of the plaquette- and 
dimer-states should emerge around j = 1 in various phys- 
ical quantities such as the excitation spectrum, etc. 

The results obtained above do not necessarily imply 
that a disordered ground state is always realized in the 
whole range of j. It is needed to study how the disordered 
phases compete with possible antiferromagnetic phases 
driven by the three-dimensional (3D) exchange couplings. 
We first recall that for j = and A = 1, the system 
should have an antiferromagnetic order, as mentioned be- 
fore. On the other hand, in the case j — > 00 ( J' — * 0) the 
spin system is reduced to the s = 1/2 massless Heisen- 
berg chain-characterized by the Tomonaga-Luttinger liq- 
uid phaseo In the parameter regime (j > 1), a different 
type of the magnetic order was observed experimentally 
for GeCu 2 04 (j ~ 6) .El Therefore, we have to carefully 
check whether the above different antiferromagnetic or- 
ders are indeed stabilized in our model. 

We first study the magnetically ordered phase in the 
region < j < 1. To this end, we compute the sus- 
ceptibility for a staggered field and the triplet excitation 
energy up to the fourth order in A for various values of 
j. To observe the second-order transition to the magneti- 
cally ordered phase, we study the spin gap at k = in the 
Brillouin zone, which should vanish at the phase transi- 
tion ppipt. By applying Pade approximants to computed 
seriesjij we obtain the phase boundaries shown in Fig. |j. 
When (j, A) = (0,0), the system is reduced to an assem- 
bly of isolated tetrahedra. With the increase of A, the 
3D network develops, enhancing the antiferromagnetic 
correlation. At last, the second-order quantum phase 
transition occurs to the magnetically ordered phase (e.g. 
the critical value is given by A c ~ 0.4 for j = 0). The 
increase of j suppresses the magnetic correlation due to 
strong frustration, and thus favors the plaquette phase. 
The critical value is estimated as j c ~ 0.9 for A = 1, 
where our generalized model is reduced to the original 
system. As seen in this figure, the boundaries determined 
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FIG. 3: Phase diagram for the s = 1/2 quantum Heisenberg 
model with pyrochlore structure. Bold line represents the 
phase boundary which separates the dimer and the plaque- 
tte phases. Solid (dashed) line indicates the phase boundary 
between the plaquette and the magnetically ordered phases, 
which is determined by the spin gap (staggered susceptibil- 

ity). 



in two distinct ways slightly differ from each other, which 
may be due to the lower-order (fourth) series expansion 
done here. Although it is desired to perform a higher- 
order calculation to determine the phase boundary more 
precisely, its essential feature is certainly given by the 
present calculation; e.g. the magnetic phase is not dom- 
inant for j = 1, but its phase boundary is rather close to 
./ I- 

We next examine another possibility _pf the antifer- 
romagnetic order observed for GeC^Oja in the region 
j > 1. For this purpose, we calculate the susceptibility 
for the corresponding staggered field up to the third order 
in A. As a result, we find that the divergent singularity 
around A = 1 is gradually suppressed, as j is decreased 
from the value (j — oo) for the isolated spin chain. This 
tendency implies that in the region j > 1, the system 
does not enter the antiferromagnetically ordered phase, 
but always stays in the dimer phase with spin gap. There- 
fore, it is seen from the above analysis that-the magnetic 
order observed for the compound GeC^Oja (j ~ 6) may 
not be simply explained in terms of the isotropic Heisen- 
berg model employed here. This in turn suggests that 
anisotropic interactions may be important to realize an 
ordered state in this compound. 

Let uS||»pw turn to the s = 1 system on a pyrochlore 
lattice.ll§ll3 Although the series-expansion calculation be- 
comes much more difficult in this case, we can still de- 
duce some instructive comments on the s = 1 pyrochlore 
lattice. In order to use series expansion techniques, we 
again start with an isolated s = 1 tetrahedron, whose 
eigenvalues are listed in Tabic g. It is seen that the 
s = 1 plaquette-singlet with S\ + S2 = S3 + S4 = 2 and 
Stotai = is the ground state for < j < 1 , whereas 
the dimer singlet state with Si + S2 = S3 + S4, = 
and iStotai = is the ground state for j > 1. In con- 
trast to the s = 1/2 model, an isolated tetrahedron in 



TABLE II: Eigenvalues of an isolated s = 1 tetrahedron: En 
(Ep) is the energy for the first (second) part in eq. (^). 



<Sl2 













1 










2 




£34 





1 


2 





1 




2 





1 






2 


Stotal 





1 


2 





1 


2 


1 2 3 


2 


1 2 


3 





12 3 4 


Eo/J 


-4 


-3 


-1 


-3 


-2 







-1 









2 


E P /J' 














-2 -1 


1 


-3 -1 2 





-3 -1 


2 


-6 


-5-3 4 



the homogeneous point (j = 1) has three-fold degen- 
erate ground states, which are composed of the above- 
mentioned singlet states together with another singlet 
state with Si + S2 = S3 + S4 = 1 and Stotai = 0, which 
may be regarded as a s = 1/2 plaquette-singlet state 
(see Table ||) . By turning on the inter-tetrahedron cou- 
pling, we observe how the above three-fold singlet states 
evolve on a pyrochlore lattice. To investigate the first- 
order quantum phase transitions among three phases, we 
estimate the ground state energy up to the fourth or- 
der in A. This expansion claims that in contrast to the 
s = 1/2 case, there exists an intermediate (s = 1/2 pla- 
quette) phase between the s — 1 plaquette phase and the 
dimer phase in small A: two phase boundaries are esti- 
mated as j c — 1 - 0.42A 3 and j c = 1 + 0.084A 3 , which are 
shown as the bold dashed lines in Fig. [|. Although we do 
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FIG. 4: Schematic phase diagram for the s = 1 antiferromag- 
netic Heisenberg model with pyrochlore structure. 

not have a definite answer to the question whether this 
intermediate state can survive as the ground state even 
for A = 1, it may be a possible candidate which com- 
petes with other states at fully frustrated point j = 1 
with A — 1. We next calculate the spin gap to see how 
stable the magnetically ordered phase for the s = 1 case 
is in comparison with the s — 1/2 case. By applying 
the Pade approximants to the third-order results Ji3 we 
deduce the phase boundary shown as the solid line in 
Fig. ^. It is seen that the area of the magnetically or- 
dered phase is more extensive than the s = 1/2 case. By 
recalling the phase diagram for the s = 1/2 case, it is 
thus expected that the magnetically ordered phase may 
be more dominant around the homogeneous point j = 1. 
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Finally, we make a brief comment on the small J' (large 
j) case. For J' = 0, the system is reduced to the s = 1 
spin chain with bond alternation, where the dimer phase 
and the Haldanc phase are separated at the critical point 
A c = 0.6.E.D Although it is difficult to estimate the phase 
boundary between these spin-gap states in the presence 
of interchain coupling, it is naively expected that the Hal- 
dane phase may disappear when j is decreased down to 
J = l- 

In conclusion, we have discussed the ground-state 
phase diagram for the s = 1/2 Heisenberg model with 
pyrochlore structure by means of the series expansion 
method. In particular, it has been found that the two dif- 
ferent spin-gap states strongly compete with each-.other 
around j = 1, where the compound Y(Sc)Mn2u may 
be located. Also, the antiferromagnetic phase has been 
shown to be extended rather,closely to the phase bound- 
ary. Concerning GeC^OjS, for which j ~ 6, it has 
turned out that the present model may not describe its 



magnetic order, suggesting that some other mechanism 
should be considered for the magnetism. For the s = 1 
system, we have not been able to deduce the definite 
conclusion on the phase diagram, but have checked that 
the magnetically ordered phase may be more dominant 
around j = 1 in comparison with the s — 1/2 case. Also, 
besides the known states such as dimer, plaquette and 
magnetically ordered states, another intermediate spin- 
gap state may also be a candidate for the ground state 
around j = 1. Since this argument has been based on the 
calculation for small A, it is desired to confirm whether 
this spin-gap state really takes part in the strong compe- 
tition around j = 1, which is now under consideration. 
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